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Abstract 

We study an one-hop device-to-device (D2D) assisted wireless caching network, where popular 
files are randomly and independently cached in the memory of end-users. Each user may obtain the 
requested files from its own memory without any transmission, or from a helper through an one- 
hop D2D transmission, or from the base station (BS). We formulate a joint D2D link scheduling and 
power allocation problem to maximize the system throughput. However, the problem is non-convex and 
obtaining an optimal solution is computationally hard. Alternatively, we decompose the problem into a 
D2D link scheduling problem and an optimal power allocation problem. To solve the two subproblems, 
we first develop a D2D link scheduling algorithm to select the largest number of D2D links satisfying 
both the signal to interference plus noise ratio (SINR) and the transmit power constraints. Then, we 
develop an optimal power allocation algorithm to maximize the minimum transmission rate of the 
scheduled D2D links. Numerical results indicate that both the number of the scheduled D2D links 
and the system throughput can be improved simultaneously with the Zipf-distribution caching scheme, 
the proposed D2D link scheduling algorithm, and the proposed optimal power allocation algorithm 
compared with the state of arts. 
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I. Introduction 

With the emergence of various applications, especially multimedia, mobile data has been 
explosively increasing in recent years. It is predicted that the mobile traffic will grow by more 
than 200 to 2000 times in the next few years according to a study from IMT 2020 [HI . The huge 
data amount pushes operators to provide high-throughput wireless access services. However, the 
current wireless access technology is already quite close to the performance limits and thus new 
communication strategies are needed to meet the increasing requirement from mobile subscribers. 

One promising approach is to use heterogeneous networks [|3, O, where one cell is divided 
into multiple small cells, i.e., microcells, picocells, and femtocells. Within each small cell, one 
low-power base station (LPBS) is equipped to serve the users in coverage. Specially, each 
LPBS is connected to the base station (BS) of the cell with a backhaul, e.g., high-speed 
fiber dH. Then, the requested files by a user are first transmitted from the BS to the LPBS 
through the backhaul and then transmitted from the LPBS to the user. The deployment of LPBS 
shortens the wireless transmission distance and results in a spatial gain, which boosts the system 
throughput dramatically compared with the conventional cellular networks. However, it may be 
very expensive to establish and maintain the LPBS as well as the backhaul. This hinders the 
application of the heterogeneous networks in practice. 

Alternatively, wireless caching recently has attracted a lot of attentions for the advantages of 
fast response and without heavily relying on the backhaul dH, dll. As observed that some files, 
e.g., video clips, usually remain popular in a certain period of time, say one day, one week, 
even one month. Meanwhile, with the fast development of integrated circuit (IC) technologies, 
the price of storage memory drops quickly. Then, the storage capacity in devices can be utilized 
for wireless caching. Specifically, some popular files can be cached into the memory of users 
during the off-load time, say middle night, such that users may obtain the requested files from 
the devices rather than the BS in the peak-load time, say daytime. Notably, a popular file may 
be requested by multiple users. Thus, more user requests can be accommodated with the help of 
wireless caching during the peak-load time. This increases the system throughput significantly. 
Besides, high-speed backhauls are not required for the wireless caching since the cached files 
can be transmitted through wireless channels during the off-load time. 

In this paper, we focus on the wireless caching networks with one-hop device-to-device (D2D) 
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communication among users, since D2D eommunication has been shown to be a promising 
candidate to improve the system throughput |I71. Then, some popular files are eaehed into the 
memory of users and eaeh user may obtain the requested files from its own memory without 
any transmission, or from a helper through an one-hop D2D transmission, or from the BS. 
We formulate a joint D2D link seheduling and power allocation problem to maximize the 
system throughput. However, the problem is non-convex and obtaining an optimal solution is 
computationally hard. Alternatively, we seek to obtain a suboptimal solution with reasonable 
complexity. Briefly, we intend to sehedule the D2D links with strong eommunieation ehannels 
and weak interference channels and allocate the power to the scheduled D2D links fairly. Thus, 
we deeompose the joint optimization problem into a D2D link seheduling problem and an 
optimal power alloeation problem. To solve the two subproblems, we first develop a D2D link 
seheduling algorithm to seleet the largest number of D2D links satisfying both the signal to 
interference plus noise ratio (SINK) and the transmit power constraints. Then, we develop an 
optimal power alloeation algorithm to maximize the minimum transmission rate of the seheduled 
D2D links. Numerieal results indieate that both the number of the seheduled D2D links and the 
system throughput ean be improved simultaneously with the Zipf-distribution caching scheme, the 
proposed D2D link seheduling algorithm, and the proposed optimal power allocation algorithm 
compared with the state of arts. 

Related literature: Wireless caching has recently been studied from theoretical perspectives in 
ll8l- [IT3ll . Speeifically, |[8l and dU study a typieal server-user model, where a file server transmits 
data bits on a shared link (broadeasting channel) to satisfy the request of eaeh user. The objeetive 
is to minimize the transmission rate on the shared link. Then, |[8l and @ propose centralized and 
decentralized eaehing sehemes to exploit both local and global caching gains, and thus achieve 
multiplieative peak rate reduetion on the shared link. ifTOl proposes an online coded eaehing 
seheme and aehieves the order-optimal long-term performanee on the shared link. Besides the 
single-layer caching in [|^- lfT0l . ifTTI proposes an order-optimal file plaeement and delivery 
seheme for a two-layer wireless caching network. [fT^ generalizes the result in [fTTll and studies 
the multi-requests wireless eaehing networks, where eaeh user requests more than one files. 
Then, lIT^ proposes a eaehing seheme based on multiple groupeast index eoding and aehieves 
the order-optimal performance. ifTSll extends the result in [fT^ and provides the complete order- 
optimal eharacterization of the transmission rate on the shard link. 
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Considering the popularity diversity of files, IfT^ - lfTTII study nonuniform eoded eaehing net¬ 
works, where files are assumed to have different popularity. [|T4l eonsiders a wireless eaehing 
network with one helper and multiple users, and develops an order-optimal eoded eaehing 
seheme. IfTSl generalizes the result in lfT4l to multiple helpers and users, and optimizes the 
trade-off among the peak rate on the shared link, the memory size in helpers, and the aeeess 
eost of users. [fT^ divides the popularity of files into several diserete levels and derives an 
information-theoretie outer bound for the nonuniform network. Different from ffT^ . where the 
peak rate on the shared link is eonsidered, lUTll optimizes the long-term performanee. That is, 
[fTTll eonsiders the average rate on the shared link and develops simple order-optimal sehemes. 
Furthermore, [|T^ - lf2^ study more praetieal seenarios and design effeetive eaehing sehemes 
subjeet to heterogenous eaehe sizes, delivery delay eonstraint, seeurity problem, prieing problem, 
and streaming sehedule problem. 

To further improve the system performanee, D2D eommunieation is proposed for wireless 
eaehing [|^ - lf2^ . Il23l eonsiders the D2D eommunieation in wireless eaehing networks from 
the perspeetive of information theory and proposes a deterministie eaehing seheme and a random 
eaehing seheme, both of whieh may aehieve the information theoretie outer bound within 
a eonstant multiplieative faetor. Then, [|2^ and ll25l provide the basie prineiple and system 
performanee of a wireless eaehing network with D2D eommunieation, and show that the gain 
from the unieast D2D eommunieation is eomparable to the gain from the eoded BS multieast. 
[|2^ proposes a novel arehiteeture to inerease the system throughput and obtains the optimal 
eollaboration distanee of D2D eommunieation. Different from [|^ - If26l . where one-hop D2D 
eommunieation is allowed, [f27ll and [|2^ eonsider multi-hop D2D wireless eaehing networks, 
where a requested file ean be direetly obtained through an one-hop D2D transmission or indireetly 
obtained after multi-hop D2D transmissions. More speeifieally, ifZTl and [|2^ study the throughput 
sealing law and propose a deeentralized eaehing seheme and a multi-hop D2D transmission 
seheme. With the sehemes, the optimal throughput sealing law is aehieved and outperforms the 
sealing law in one-hop D2D eommunieation networks. 

Contributions: We adopt the arehiteeture similar to [l26ll . where eaeh user may obtain the 
requested files from its own memory without any transmission, or from a helper through an one- 
hop D2D transmission, or from the BS. However, our work has three main different aspeets from 
the previous work in [|2^ . Firstly, we eonsider a general eonventional model, where any two elose 
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users are allowed to establish a D2D link, provided that the requested file of one user is eaehed 
in the memory of the other user. This is different from the eluster-based model in [l26ll . where a 
eell is divided into multiple non-overlapping elusters and only the users in the same eluster are 
allowed to establish a D2D link. Thus, our proposed algorithm has a relaxed eonstraint on the 
loeation of the users in the eell and ean potentially ereate more D2D eommunieations among 
users. Seeondly, we manage the mutual interferenee among different D2D links by effieiently 
seheduling the D2D links and fairly alloeating the transmit power. However, Il26l allows at most 
one D2D link to work in eaeh eluster to avoid strong interferenee. Thirdly, we eonsider the 
max-min fairness of the seheduled D2D links to aehieve perfeet fairness and study the optimal 
power alloeation at D2D transmitters. This is different from ll26l . where the fairness is ignored. 
As a result, the three differenees lead to a eompletely different problem formulation and solution 
eompared with that in If26l . To further elarity, our eontributions are listed as follows. 

• We eonsider the eonventional model, where any two elose users are allowed to establish a 
D2D link provided that the requested file of one user is eaehed in the memory of the other 
user. This ereates more D2D eommunieations among users eompared with the eluster-based 
model in [|2^ . 

• We manage the mutual interferenee among different D2D links by effieiently seheduling 
the D2D links and fairly alloeating the power alloeation instead of dividing a eell into 
non-overlapping elusters and allowing at most one D2D link in eaeh eluster in [|^ . 

. We formulate a joint D2D link seheduling and power alloeation problem to maximize 
the system throughput. Due to the non-eonvexity of the joint optimization problem, we 
deeompose it into a D2D link seheduling problem and an optimal power alloeation problem. 
Speeifieally, we intend to first maximize the number of the seheduled D2D links by solving 
the D2D link seheduling problem and then maximize the minimum transmission rate of the 
seheduled D2D links by solving the optimal power alloeation problem. 

. Numerieal results indieate that both the number of the seheduled D2D links and the system 
throughput ean be improved simultaneously with the Zipf-distribution eaehing seheme, 
the proposed D2D link seheduling algorithm, and the proposed optimal power alloeation 
algorithm eompared with the eaehing sehemes and D2D link seheduling algorithms in ll25l 
and [[26ll . 
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II. System Model 

We consider a cellular network with one BS and K users, where the BS is located at the 
center of a cell and the users are randomly distributed in the cell as shown in Fig. [B The BS 
has a file sever and stores N files with equal sizeQ, whose popularity probability follows Zipf 
distribution ll29l . i.e., ^ 

/.= ( 1 ) 

where t] is the file index, and 7 ^ is the file request coefficient and controls the popularity 
distribution of files. Namely, a large 7 ^ means that the first few files dominate the requests from 
users. 



Base Station 
(BS) 

User 


Figure 1. System model, where a user can obtain the requested file from its own memory (dash arrow), or from a helper 
through an one-hop D2D transmission (dash-dot arrow), or from the BS (solid arrow). 


The file distribution by wireless caching consists of a placement phase and multiple delivery 
phases. In the placement phase, each user randomly and independently caches one out of N 


‘This can be justified by the fact that big files are broken into small fragments with the same length. 
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files in its memory, according to the Zipf-distribution with a file caching coefficient 7c 111 In 
each delivery phase, each user requests one file. If the requested file of a user can be found 
in its own memory, the user accesses the file without any transmission. On the other hand, if 
the user cannot find the requested file in its own memory, it can request the file either from a 
helper through an one-hop D2D transmission or from the BS. Here, any user in the cell can be 
a helper of a requesting user if two conditions are satisfied. One is that the two users are close 
in the cell. This is, the requesting user is in the help region, which is defined by a help distance 
r around the helper. This may guarantee small pathloss and small transmit power of a D2D 
link on average. The other one is that the requested file is cached in the helper’s memory. It 
should be noted that, a user may be served by more than one helpers. Then, the helper that can 
provide the highest transmission rate or SINK will be chosen as a transmitter in a D2D link. For 
clarity, the self-access, the one-hop D2D transmission, and the BS transmission are illustrated as 
dash arrow, dash-dot arrow, and solid arrow in Fig. [H respectively, and the direction of a arrow 
denotes the transmission direction of a file. 

We assume that all the D2D links share one WIFI channel and each D2D user works in 
half-duplex mode and receives or transmits wireless data in one time slot. We also assume that 
each user transmits data in the unicast mode and then cannot transmit data for multiple users 
simultaneously. This can be justified by two facts. One is that the asynchronous nature of the 
requested files, e.g., video clips. This means that the probability that different users request the 
same video clip from a common helper at the same time is negligible llTdll . The other one is 
that the mobile nature of users. This makes it very hard to adopt the coded caching schemes 
which exploit the multicast opportunities between one transmitter and multiple receivers in |[8l 

^We consider one file to simplify the presentation, although it is straightforward to generalize to multiple files with a Zipf 
distribution by caching each file independently. Note that the Zipf-distribution caching is not optimal for caching one file or 
multiple files. We use the Zipf distribution to cache files for four reasons. The first one is that, we focus on the algorithm 
design in the delivery phase and thus the proposed algorithms in this paper can be used for any caching scheme. The second 
one is that, the optimal caching distribution of the general model is quite hard to obtain (if it is possible). This is because the 
optimal caching scheme of one user is related to the number of its neighbors, which is defined as the users in its help region. 
In (21 and 1251 . a cell is divided into multiple non-overlapping clusters and the users in each cluster are assumed to be the 
neighbors of each other. Then, the numbers of neighbors for all the users in the same cluster are identical. This simplifies the 
mathematical derivation and makes it possible to obtain the optimal caching distribution for the cluster model. In our general 
model, the numbers of neighbors in terms of different users are quite different, especially for the users in the center of a cell and 
the users at the edge of a cell. This makes it quite hard to analyze the optimal caching distribution mathematically for a general 
model. Thus, a further study on the optimal caching scheme is beyond the scope of this paper. The third one is that, the Zipf 
distribution can well model the popularity of files 1291 . Then, we use the Zipf distribution for wireless caching as a heuristic 
choice. The forth one is that, from the numerical results, the system performance (both the number of the scheduled D2D links 
and the system throughput) with the Zipf-distribution caching scheme and the proposed algorithms in this paper outperforms 
the system performance with the optimal caching distribution for the cluster model and the algorithms in (in and (25 1 . 


and O . In fact, the request of a video elip in most praetieal systems is implemented by a point 
to point (a transmitter to a reeeiver) transmission with a dedieated eonneetion [|2^ . 

In the following, we will evaluate the SINK at D2D reeeivers. Suppose that ks users in set Ss 
ean aeeess their requested files from their own memory, and /cdb users in set iSdb can obtain the 
requested files either from helpers through one-hop D2D transmissions or from the BS, and 
users in set iSb can only obtain the requested files from the BS. Define the link through whieh 
a user in iSdb can obtain the requested file by D2D transmission as a potential D2D link and 
denote Im as the potential D2D link with reeeiver m. Without loss of generality, if we assume 
^DB = {1,2, • • • , /cdb}, the potential D2D link set ean be denoted as £db = hr " > hm}- 
Suppose that all the users in iSdb are seheduled to work simultaneously and that the ehannel 
gain between the transmitter of the D2D link {rn G iSdb) to the reeeiver of the D2D link 
(n G iSdb) is g{m,n). Then, the SINK at the user m (D2D reeeiver) is 

Pmg{m,m) 

J:n,S..,n^mPn9in,m) + Nj 

where pm is the transmit power of the D2D link Im and it is subjeet to the transmit power 
eonstraint 0 < pm < and Nm is the power of the additive white Gaussian noise (AWGN) 
at the reeeiver m. Here, we eonsider quality of experience (QoE) guaranteed D2D transmissions. 
Thus, the SINK at eaeh D2D reeeiver is required to be greater than or equal to the minimum 
aeeeptable SINK, i.e., Vm > Vm, whieh is related to the modulated and eneoding seheme at eaeh 
D2D transmitter. 


III. Problem Formulation and Analysis 
A. Problem Formulation 

From the previous seetion, there is a ehanee that not all the potential D2D links in £db can 
be seheduled simultaneously for two reasons. One reason is that the potential D2D links in 
Fdb may share the same users. That is, one user may be a D2D transmitter in one potential 
D2D link and a D2D reeeiver in another potential D2D link. Then, the two links eannot work 
simultaneously. The other reason is that, some potential D2D links eannot be satisfied with the 
minimum aeeeptable SINRs. Thus, we intend to develop a joint D2D link seheduling and power 



allocation algorithm to maximize the system throughput. Formally, we have 


(Po): 


max W log(ti„ + 1) 

Pm :m.60D ^—' 


S.t. 


iSd c iS] 


DB; 


iSdt n iSd — 0, 

Vm>Vm,^me 5 d , 

0 < Pm < Pm^^V m e 5 d, 


(3) 

(4) 

(5) 

( 6 ) 


where the seheduled D2D link set iSd and the eorresponding transmit power p^, V m 6 iSd, are 
the variables to be optimized, iSdt is the D2D transmitter set eorresponding to the D2D reeeivers 
in iSd, and the eonstraint (4) guarantees that different seheduled D2D links do not share the same 
D2D users. 

Unfortunately, the objeetive funetion of (Pq) is non-eonvex. To obtain the joint optimal 
D2D link seheduling and power alloeation solution, the eomplexity of the exhaustive seareh 
is quite high 0. Alternatively, we seek to obtain a suboptimal solution of (Pq) with reasonable 
eomplexity. Briefly, we intend to sehedule the D2D links with strong eommunieation ehannels 
and weak interferenee ehannels and alloeate the power to the seheduled D2D links fairly. Thus, 
we deeompose (Pq) into two subproblem^ i.e., a D2D link seheduling problem and an optimal 
power alloeation problem. For the first subproblem, we intend to optimize iSd, to maximize the 
number of D2D links satisfying both the SINK and the transmit power eonstraints. For the seeond 
subproblem, we intend to optimize pm, V m G iSd, to maximize the minimum transmission rate 
of the seheduled D2D links. In what follows, we will formulate the two subproblems. 

1) D2D Link Scheduling Problem: To maximize the number of the seheduled D2D links 
satisfying both the SINK and the transmit power eonstraints, the D2D link seheduling problem 


, / _, max \ 

^The complexity of the exhaustive search is O 1, where is the maximum transmit power of the D2D 

link Im, e is, the maximum tolerance error of the exhaustive search scheme, and fcoB is the number of users that can obtain the 
requested files either from helpers through one-hop D2D transmissions or from the BS. Suppose that /cdb = 20, = 

pinax _ 2 qq V m 7 ^ n, and e = 0.01. Then, the complexity of the exhaustive search is extremely large, i.e., O (10000^°). 

®Note that the decomposition is not optimal and finding the optimal decomposition will be our future work (if it is feasible). 
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can be written as 


(Pi) : max IiSdI 
<SdC<Sdb 

s.t. (4), (6), 

t'm > Wm, V m G 5 d, (7) 

where |iSd| is the cardinality of iSd and Vm = It is clear that the solution of (Pi) may maximize 
the number of the scheduled D2D links. However, scheduling the largest number of D2D links 
may lead to small system throughput since more D2D links may cause more mutual interference. 
To deal with this issue, we introduce a scheduling coefficient Cg (dB) and replace Vm = with 
Vm = niaxlfim, Cs}- Then, if c* < min Vm, we have Vm = Then, (Pi) maximizes the number 

meSoB 

of the scheduled D2D links. As we increase Cg, the SINK constraint (7) becomes strict, only the 
D2D links with strong communication channels and weak interference channels are scheduled. 
This potentially increases the system throughput. Thus, we can numerically choose the value of 
Cs to obtain an acceptable number of the scheduled D2D links as well as the system throughput. 

Theorem 1: Each transmit power pn, ^ n e iSd, increases if any SINK Vm, rn G iSd, in (Pi) 
increases. Meanwhile, the value of |iSd| remains constant or decreases if any SINK Vm, rn G iSd, 
in (Pi) increases. 

Proof: The proof is provided in Appendix A. ■ 

From Theorem 1, the value of |iSd| may decrease if any SINK Vm in iSd increases. Then, 
problem (Pi) is equivalent to 


(P 2 ) : max |5 d| 

‘SdCSdb 

S.t. (4), (6), 

Vm = Vm^ rn e Sji. (8) 

By solving problem (P 2 ), we may obtain the optimal D2D receiver set and the correspond¬ 
ing D2D link set Without loss of generality, we assume = {1, 2, • • ■ , /cd} and = 
{^ 1 ,^ 2 , • • • ,4d}- It is clear that problem (P 2 ) outputs the SINK vector V = [fii,h 2 , • • • 
at the scheduled D2D receivers. Accordingly, we assume the power allocation vector is P = 
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2) Optimal Power Allocation Problem: From Theorem 1 again, the value of |iSd| may remain 
eonstant if any SINK Vm, m G S^, inereases. Thus, we may improve the minimum transmission 
rate or the minimum SINK of the seheduled D2D links without eompromising |iSq|, provided 
that the transmit power constraints of the scheduled D2D links are not violated, i.e.. 


P 3 ) : max min \og{vm + 1) 

Pm'-mCLS^l<m.<kD 


s.t. Vm > m e S^, 

(9) 

Pm<Pm< V m G iSp. 

( 10 ) 


B. Analysis of Subproblem {P 2 ) 

To solve problem (P 2 ), we shall develop a D2D link scheduling algorithm to select the 
largest number of D2D links satisfying all the constraints in problem (P 2 ). In fact, problem 
(P 2 ) is an admission control problem (i.e., link scheduling problem) [l^ - lf35l . However, it is 
different from a regular admission control problem. In the regular admission control problem, 
one transmitter has one dedicated receiver. In our system model, there is a chance that one user 
is the transmitter in one potential D2D link and the receiver of another potential D2D link. Since 
the scheduled D2D links cannot share the same D2D users, i.e., the constraint (4), problem (P 2 ) 
is more complicate than the regular admission control problem and the algorithm for the regular 
admission control problem cannot be directly used to solve problem (P 2 ). Thus, we need to 
develop a scheduling algorithm to solve problem (P 2 ). 

Denote the events that any two potential D2D links in a potential D2D link set do not share 
the same users and that all the minimum SINK constraints in a potential D2D link set can be 
satisfied as Ci and C 2 . Accordingly, Ci and C 2 denote the events that Ci and C 2 cannot be 
satisfied, respectively. Then, we have the following three cases depending on whether Ci and/or 
C 2 can be satisfied in the potential D2D link set £db- 

Case I: (Ci, 02 ). In this case, all the potential D2D links in £db will be scheduled. Then, we 
have = Sub- 

Case II: (Ci, C 2 ). In this case, problem (P 2 ) reduces to a regular admission control problem. 
To solve this problem, the fewest D2D links in £db shall be removed until C 2 can be satisfied. 

Case III: (Ci). In this case, some D2D links share the same users and are not allowed to work 
simultaneously. To solve this problem, we shall first divide the potential D2D links in £db into 
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different groups (or subsets), such that any two potential D2D links in one group do not share 
the same users. After that, each group becomes case I or case 11. 

For generality, we will consider Case III for problem (P 2 ) in the rest of this paper. By analyzing 
problem (P 2 ), we have 

Theorem 2: Problem (P 2 ) is NP-hard. 

Proof: The proof is provided in Appendix B. ■ 

Considering the highly computational complexity of obtaining the optimal solution for a NP- 
hard problem, we will develop an efficient algorithm to obtain a suboptimal solution of problem 
(P 2 ) in what follows. 

C. Analysis of Subproblem {Pfj 

Once the optimal D2D set is obtained from problem (P 2 ), problem (P 3 ) is a feasible 
max-min optimization problem with transmit power of D2D links as variables. Then, we will 
first analyze the property of the optimal power allocation of problem (P 3 ) and then develop an 
optimal power allocation algorithm to achieve the optimal performance. 

IV. D2D LINK Scheduling 

In this section, we aim to maximize the number of the scheduled D2D links. However, some 
of the potential D2D links may be dependent. That is, different potential D2D links in £db 
share the same users. This is the main difference between problem (P 2 ) and a regular admission 
control problem, where any two different potential D2D links do not share the same users. 
Thus, we shall divide the potential D2D links in £db into different groups, such that different 
potential D2D links in each group are independent. In this way, problem (P 2 ) is decomposed 
into several subproblems corresponding to different groups. After maximizing the number of the 
potential D2D links that satisfying both Ci and C 2 for each subproblem (group), we can obtain 
a suboptimal solution of problem (P 2 ). 

In what follows, we will first divide the potential D2D links in £db into different groups 
such that the potential D2D links in each group are independent and satisfy Ci. Then, we will 
maximize the number of the potential D2D links satisfying C 2 in each group. Specifically, we 
develop a centralized power control (CPC) algorithm to check the state of each group, i.e., 
to decide whether the potential D2D links in each group satisfy C 2 or not. If the potential 
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D2D links in one group do not satisfy C 2 , we further develop a removal algorithm to remove 
the fewest potential D2D links from this group. After applying the CPC algorithm and/or the 
removal algorithm into eaeh group, the largest number of the potential D2D links satisfying 
both Cl and C 2 in eaeh group is obtained. Finally, the potential D2D links in the group with 
the largest number of potential D2D links will be seheduled and other potential D2D links will 
not be allowed. 

A. Potential D2D Links Division 

To maximize the number of the seheduled D2D links, we shall add as many potential D2D 
links as possible into eaeh group. This leads to the fewest groups. Thus, the division of the 
potential D2D links ean be transferred to an edge eoloring problem in graph theory, where the 
edges in a graph are eolored with the fewest eolors while guaranteeing that any two edges sharing 
the same vertex are eolored differently. 

Denote the eonneetion among the users involved the potential D2D links in £db as a graph 
G{0,8), where O is the vertex set denoting the involved users in £db, 8 is the edge set 
denoting the D2D links in £db- Then, we may eolor the edges in 8 with x{G) eolors, whieh 
is the edge ehromatie number [l30ll . Sinee the edge eoloring problem is NP-hard for a general 
graph, it is non-trivial to obtain the optimal solution. Alternatively, we develop an iterative 
edge-coloring algorithm as shown in Algorithm [H which colors the edges in 8 with Nq colors 
{Nc > x{G)). More formally, denote 8 = 8iVJ ...VJ 8 nc^ where 8i {1 <i < Nq) represents the 
edge set with the Ah color after edge coloring. We assume that there are ki edges in 8i, which 
corresponds to potential D2D links. If we denote the receiver set of the ki potential D2D links 
as Si = {A, * 2 , • • ■ 5 the ki potential D2D link set can be denoted as li 2 ,- " ^ }• 

In the Ah iteration of Algorithm [H the maximum matching 0 is obtained by solving a linear 


^The maximum matching is the maximum edge set, where any two edges do not share a common vertex. 
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programming problem OTI 



S.t. ^ Xe < 1,V O e (9, 


( 11 ) 


0 < Xe < 1, V e G 


( 12 ) 


where e ~ o means that the edge e is ineident on the vertex o, and rounding the solutions Xg 
(e G £^) to 0 or 1 in order to map the solutions to integers. That is, Xg = 0 if 0 < Xg < 0.5 and 
Xg = 1 if 0.5 < Xg < 1. In this way, the edges Xg = 1 are eolored with the zth eolor and put 
into the set = {/(e) : Xg = 1, e G £^}, where /(e) denotes the D2D link eorresponding to the 
edge e. This proeedure is terminated until all the edges are eolored. Note that the eomputational 
eomplexity of Algorithm [T] is dominated by solving the linear programming problem (P 4 ) with 
eomputational eomplexity 0 (/^db) IEH- Thus, the overall eomputational eomplexity of Algorithm 


His 0 (/cdb)- 


Algorithm 1 Iterative Edge-eoloring Algorithm. 

Initialization 
1 : G{0,8y, i= 0 ; 

Iterative: 

2 : while 8 do 
3: i = -/ + 1, = 0; 

4 : Solve (P 4 ) and obtain Xg, V e G 8, and round Xg to 0 or 1; 

5: Color all the edges Xg = 1 with the ith eolor; 

6 : Ci = {l{e) ■. Xe = l,e e 8}, 8 = 8\{e : Xg = 1}, ; 

7: end while 


B. CPC Algorithm 

After dividing the potential D2D links into Ci < i < Nq) with the iterative edge-eoloring 
algorithm, Ci is satisfied among the potential D2D links in eaeh In what follows, we will 
develop an algorithm to eheek the state of the potential D2D links in £j, i.e., to deeide whether 
the potential D2D links in Ci satisfy C 2 or not. 

If all the potential D2D links in Ci satisfy C 2 , there exists a power alloeation to satisfy all 
the potential D2D links in Ci with the minimum SINR eonstraints Vj = Vi^^ ..., Vi^Y■ 
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Suppose that the power allocation is Pj = Pi^^Y- Then, we have 


//,(V,)P, = N„ 


( 13 ) 


where 


T 




-g{i2,il) 5'(*2,*2)M2 


g{iii '^ki) 

g{i2,iki) 


(14) 


g{iki)ii) giiki)i2) ■■■ g{'^ki,iki)/vif.. 


and 



T 


(15) 


Then, the power allocation Pj can be derived as 




(16) 


If the potential D2D links in £j satisfy C 2 , we have 0 ^ P* ^ P™, where P™'™ = 
^ y ''' yPT^Yy 0 is a fcj X 1 vcctor with all zero elements, and is an element- 
wise operator. 

Consequently, we have the CPC algorithm as 

• the potential D2D links in Ci satisfy C 2 if 0 ^ Pj ^ pmax 

• the potential D2D links in £j do not satisfy C 2 if 0 ^ Pj ^ pmax holds. 

In our algorithm, the BS first collects the channel gains among different potential D2D links 
and the minimum SINK constraint of each potential D2D link, and then calculates the required 
power allocation to satisfy the potential D2D links in Ci with their minimum SINK constraints. 
Thus, we name this algorithm as centralized power control (CPC) algorithm. Essentially, the main 
ideas of the CPC algorithm and the distributed constrain power control (DCPC) algorithm in lf32l 
are similar: check whether multiple transceivers can be supported simultaneously by checking 
whether the transmit power constraints are violated to achieve the minimum SINK constraints. 
However, the CPC algorithm is more efficient than the DCPC algorithm to implement. More 
specifically, to implement the DCPC algorithm, each D2D transmitter in Ci iteratively adapts its 
transmit power to achieve the minimum SINK constraint until all the transmit power converges. 
If all the minimum SINK constraints are satisfied after transmit power converges, the potential 
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D2D links in £j satisfy C 2 . If some minimum SINK constraints are not satisfied, the potential 
D2D links in do not satisfy C 2 . Then, a removal algorithm will be developed to remove 
some potential D2D links from In the CPC algorithm, the BS first collects all the necessary 
information and then calculates the power allocation to achieve the minimum SINK constraints. 
If the calculated power allocation can be satisfied at all the D2D transmitters in the potential 
D2D links in satisfy C 2 . Otherwise, a removal algorithm will be developed. 


C. Removal Algorithm 


From the CPC algorithm, there exists at least one calculated transmit power pi^ in £* satisfying 
Pim > pST Pim < 0 if the potential D2D links in do not satisfy C 2 . This is because the 
mutual interference among the potential D2D links in £* is too strong. Then, we shall remove 
some potential D2D links from to enable the remaining ones to satisfy C 2 . Our approach is 
to remove the potential D2D link which is likely to cause the strongest interference to others or 
receive the strongest interference from others each time until the remaining potential D2D links 
satisfy C 2 . 

Specifically, to satisfy the minimum SINK constraint (im G Si), the transmitter of the 
potential D2D link has to set the transmit power no less than and causes no less than 

in) (in ^ Si, in ^ im) interference to another potential D2D link Since a potential 
D2D link with a smaller minimum SINK constraint and a larger maximum transmit power can 
tolerate more interference from other potential D2D links, we define the relative interference 
from P^ to as Ir{P,^,P,P) = in)- Then, the summation of relative interference 

generated by P^ is larger than 


n=ki 

~ ^ ^ Pr(Pml(in) ~ 

n=l,n7^m 


Ni Vi 

Lm 


n=ki 

n=l,n^m 


Pi 


79(im-i ir, 


(17) 


Similarly, to satisfy the minimum SINK constraint Vi^ (in € Si, in f im), the transmitter of the 
potential D2D link p^ has to set the transmit power no less than . Then, the summation 

of relative interference to the potential D2D link p^ is larger than 




n=ki 

z 

n=l,n^m 


Pr(Pn, (im) ~ 


Vi. 


Pt. 


n=ki __ 

NinVir 

,max g(in,ir. 


~9(in, ir. 


( 18 ) 
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Thus, the potential D2D link , where 

W = arg max max{ai^,(19) 

l<m<ki 

is likely to eause the strongest interferenee to others or reeeive the strongest interferenee from 
others and will be removed from 

It should be noted that, the proposed removal algorithm is different from the DCPC-based 
removal algorithm in [|33. Speeifieally, our proposed algorithm utilizes the information of the 
maximum transmit power at eaeh transmitter while the DCPC-based removal algorithm in 
utilizes the eonverged transmit power of the DCPC algorithm, although the two algorithms share 
other two kinds of information, i.e., the minimum SINK eonstraints and the interferenee ehannel 
gains among different potential D2D links. Besides, our proposed removal algorithm depends 
on the relative interferenee among different transmissions whereas the DCPC-based removal 
algorithm in [|^ is developed by measuring the absolute interferenee after the DCPC algorithm 
eonverges. Thus, the ealeulated interferenee in our proposed algorithm is more aeeurate than 
that in [|3^ and our proposed removal algorithm is more flexible than the DCPC-based removal 
algorithm in [[32. In faet, from the numerieal results in Seetion VI, our proposed removal 
algorithm outperforms the DCPC-based removal algorithm. 

D. Scheduling Algorithm and Complexity Analysis 

After applying the CPC algorithm and/or removal algorithm into eaeh potential D2D link set 
(1 < / < /Vc), the potential D2D links in eaeh Ci satisfy both Ci and C 2 . Then, the potential 
D2D links in = arg max |£j| with the largest number of the potential D2D links will be 

l<i<Nc 

seheduled. To summarize, we illustrate the detail D2D link seheduling algorithm in Algorithm 

m 

Note that Algorithm|2]consists of AlgorithmlH the CPC algorithm, and the removal algorithm. 
The eomputational eomplexity of Algorithm [T] is 0(/sdb)- The eomputational eomplexity of the 
CPC algorithm is dominated by ealeulating (fT^ . i.e., Pj = In the worst ease, 

is a /cdb x ^db matrix and Nj is a /cdb x 1 veetor. The eomputational eomplexity of 
the inverse of a k^B x /cdb rnatrix Hi(Vi) is 0 (/cdb^^) *^he eomputational eomplexity of the 
multiplieation of a /cdb x /cdb niatrix a /cdb x 1 veetor Nj is (/(/cdb) lEH- Then, 

the eomputational eomplexity of the CPC algorithm is 0 (/cd|^^) + 0{k‘^^) = 0{k‘^^'^^). Besides, 


18 


the computational complexity of the removal algorithm is dominated by arithmetic and thus is 
O(fcDB)- Consequently, the computational complexity of Algorithm |2] is 0(4 b) + 0(4?“) + 

0(k^„) = 0 ( 4 h ). 


Algorithm 2 D2D Link Scheduling Algorithm. 

Initialization 

1 : 0(0, S); 

Iterative: 

2: Execute Algorithm [U to color S with Nq colors, i.e., £ = £iU ... U £nc '■> 

3: for i = 1 to i = Nc do 

4 : Adopt CPC algorithm to check whether the potential D2D links in corresponding to 

the edge set £i satisfy €2', 

5 : if the potential D2D links in Ci do not satisfy C 2 then 

6 : while the potential D2D links in Ci do not satisfy C 2 do 

7 : Remove the potential D2D link im* according to (fT9l) : 

8 : Update Ci = Ci\li ^„; 

9: end while 

10: end if 

11: end for 

12: BS schedules the potential D2D links in set Ci* = arg max \Ci\. 

l<i<Nc 


V. Optimization of Power Allocation 

In the previous section, we have solved problem (P 2 ) and obtained the scheduled D2D link 
set Ci* = Cq = {h, h, - ■ ■ , 4d} and the D2D receiver set = {1,2, • • • , /cd}. Accordingly, we 
may obtain the output transmit power vector and the SINK vector as Pj* = P = [pi,p 2 -i - ■ ■ , Pk^ 
and Vj* = V = [vi,V 2 ,--- from problem (P 2 ). In this section, we will solve problem 

(P 3 ) to maximize the minimum transmission rate or the minimum SINK of the scheduled D2D 
links without compromising the number of the scheduled D2D links. Although problem (P 3 ) 
can be transformed to a convex problem and is solved with a convex optimization toolbox, we 
seek to solve it analytically to shed more light on the optimal power allocation. Specifically, we 
will first analyze the property of the optimal solution of problem (P 3 ). Then, we will develop a 
binary-search based power allocation to obtain the optimal power allocation. Finally, we analyze 
the computational complexity of our algorithm. 
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( 20 ) 


A. Binary-search based Power Allocation 

As mentioned above, the minimum SINK constraints V = V = [hi,i; 2 ,-- - of the 

scheduled D2D links are satisfied with the power allocation P = [pi,P 2 , •' ‘ iPkoY■ To further 
increase the SINRs in V and maximize the minimum SINK in V, i.e., problem (P 3 ), we shall 
increase the SINRs in V with small minimum SINR constraints with priority. 

Intuitively, if we assume hi < h 2 < • • • < Vk^, we shall increase vi from hi to h 2 , i.e., 
hi < fi < V 2 , and then increase vi and V 2 from V 2 to ha, i.e., h 2 < fi = f 2 < hs, and then 
increase vi, V 2 , ■ ■ ■, Vm 0- < m < /cd) from Vm to h^+i, i.e., h^ < vi = ^2 = • • • = < h^+i. 
This procedure is terminated until some maximum transmit power constraints are violated. In 
this way, we maximize the minimum SINR of the scheduled D2D links. More formally, we have 
the following Theorem. 

Theorem 3: Assume hi < h 2 < • • • < h^j, and denote = [ y, • - , h^+i, • • • , Vk^], we 

m 

have V(™)(hm) = [ h^, • - ,h,^^ ,h^+i, • • • ,Vkj,] and P(V(™)(hm)) = where 

m 

is given in (Ei) and N = [iVi, N 2 , • • • , 

. For 1 < m < fco, if P ^ P(V(”^)(h^)) A P“‘‘\ where P”"’^ = P = 
holds and P A P(V*^”^i(hm+i)) ^ optimal power allocation in 

problem (P 3 ) enables the optimal SINR vector V* = - ,vl^] to satisfy V* = 

= [v\_^ Vm+l, • • • , Vko], where Vm < v* < Vm+l- 

_m 

• For m = fco, if P A A jiglds, the optimal power allocation in problem 

(P 3 ) enables the optimal SINR vector V* = [^ 1 ,^ 2 , • • • to satisfy V* = Y^'^\v*) = 
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\ v*, ■ - , v*j , where v* > Vk^. 

ko 

Proof: The proof is provided in Appendix C. ■ 

Suppose that there exists a user m (1 < m < k^f}, whieh satisfies P A ^ 

but does not satisfy P A P(V’^™^(fim+i)) ^ From Theorem 3, the optimal power allocation 
P* in problem (P 3 ) enables V* = [ n*, • - , n* , Vm+i, • • • , Vko] and can be written as 

m 

P* =//-i(V*)N, (21) 

which means that the optimal power allocation P* can be obtained by calculating //~^(V*)N. 
However, the optimal SINK vector V* is unknown to the BS. This makes it difficult to obtain 
P* directly. Alternatively, we develop a binary-search based algorithm to approach V* and then 
obtain P*. Specifically, the optimal SINK vector is V* = = \ v*, ■ — , u* , Vm+i, • • • , 

m 

where Vm ^ v* < Vm+i- Consider that each transmit power in P* increases if any SINK in V* 
increases from Theorem 1. Then, we may approach V* by applying binary search between 
and subject to the transmit power constraints. 

Similarly, if P A A p^^x Folds, the optimal SINRs can be denoted as 

V* = ,u*J , where v* > Vk^- Meanwhile, consider that the optimal SINK 

kj) 

vector V* is upper-bounded by min - < y* < min ^ Then, we 

may approach V* by applying binary search between \vk^) and ^ min 

subject to transmit power constraints. To summarize, we illustrate the binary-search based power 

allocation in Algorithm |3l 

B. Complexity Analysis 

The computational complexity of Algorithm |3] is dominated by two loops. In what follows, 
we will analyze the computational complexity of the two loop, respectively. 

For the first loop from Line 1 to Line 6, we consider the worst case that there are ky, rounds. In 
each round, the computational complexity is dominated by the calculation of P(r('”^(i;m+i)) = 
which consists of a matrix inversion operation and a matrix multipli¬ 
cation operation. Since the computational complexity of the inversion of a /cd x matrix 
is 0{k^^'^^) and that of the multiplication of a /cd x kj) matrix 
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Algorithm 3 Binary-search based Power Allocation. 

Initialization 

m = 1, maximum tolerance error: e^; 

Iterative: 

1 : while P ^ P(V*^™')(hm+i)) ^ P™ do 
2 : m = m + 1; 

3: if m > then 

4: Break; 

5: end if 

6: end while 

7: if m < fco then 

8- '^rain '^mt ^max ^mid ^maxj 

9: else 

10: 'ymin ^max min jj , riniid Tlmax? 

11 : end if 

12 : while P ^ P(V*^™^(i;niid)) ^ P™^’^ does not hold or Wmax — do 

13: riniid = 

14: if P ^ P(V(”*)(r;mid)) ^ does not hold then 

15: Tljiiax ^midj 

16 : else 

17: rijiiin Tlmidj 

18 : end if 

19: end while 

20 : return P* = P(V(”^)(nniid)); 


and a /cd X 1 vector N is a 0{kl^) llJTl . we have the computational complexity to derive 

P(V(”^)(hm+i)) is -I- 0{k‘l^) = Then, the computational complexity of the 

first loop is 0{k^^^^). 

For the second loop from Line 12 to Line 19, the round of the binary search is O (log0) [f38ll . 

min ^ .j - -—Vkr. 

m^S* ^Tn D 

where 0 = —^- - -. In each round, the computational complexity is dominated by 

the calculation of P(V*^™)(r;mid)) and thus is 0{k‘^^'^^). Then, the computational complexity of 
the second loop is log 0). 

Consequently, the computational complexity of Algorithm |3] is 0{k^^'^^) + log 0 ). 

VI. Numerical Results 

In this section, we will show the performance of our proposed algorithms from numerical 
sides. To show the advantages of the proposed algorithms, we will also give the performance 
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Figure 2. Simulation network model, where the empty circle denotes the BS and a solid circle denotes a user. 



r 


Figure 3. The number of users that can obtain the requested files from their own memory or from helpers through D2D 
transmissions, i.e., ks + fco, with different help distances r and file caching coefficients 7c, where wt = 0 dB, Cs = 0 dB, 
7 r = 0.6, and K — 100. 


comparison with the existing similar algorithms in lf25l . If26l . ll^ . and ll34l . 

In our simulation, we adopt the pathloss ehannel as y = -^x + Uq, where y is the received 
signal at the reeeiver, x is the transmit signal at the transmitter, d is the distanee between two 
transceivers, A is the wavelength of the carrier, and uq is the AWGN with power spectral density 
Nq = —170 dBm/Hz at the reeeiver. We assume that all the maximum transmit power eonstraints 
are the same, i.e., = 20 dBm, W 1 < m < K, and that the minimum aeeeptable SINRs 

at all the scheduled D2D receivers are the same, i.e., Vm = vt, V m G S^, that the number 
of overall files is iV = 1000, and the bandwidth for D2D communieation is 1 MHz. For fair 
eomparison with [l26]l . we assume a square eell with eaeh side 1 km as shown in Fig. [2l 

Fig. [3] provides the number of users that ean obtain the requested file from their own memory 
or from helpers through D2D transmissions, i.e., ks + kj), with different help distances r and 
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K 

Figure 4. The Optimal help distance versus different number of users in the cell, where vt = 0 dB, Cs = 0 dB, jr = 0.6. 


file caching coefficients 7 c. From this figure, the ks + ku first increases and then decreases 
as r grows. In fact, although ks is only determined by the caching coefficient, kj) is directly 
affected by two factors. One is the number of the potential D2D links ^db- The other one is the 
interference among the scheduled D2D links. Then, ks + k^ is more sensitive to the variation 
of ki) compared with ks- As r grows, the chance that one user finds the requested file from 
other users memory increases. This increases /cdb, and ks + ki). As r continues to grow, the 
average distance between D2D transceivers increases. This leads to higher transmit power at 
each D2D transmitter. Thus, the mutual interference among the scheduled D2D links increases, 
and decreases ks + kj). Besides, we observe that the optimal r and 7 c in terms of the largest 
ks + ku is about 0.33 km and 1.5, respectively. 

Fig. m gives the optimal r with different number of users K. For comparison, we also provide 
the optimal r in [[26]l . It is clear that the optimal r decreases as K grows. For small K, the 
number of the potential D2D links /cdb limits the number of users that can obtain the requested 
file from their own memory or from helpers through D2D transmissions, i.e., fcs+^o- To maximize 
ks + ku, the optimal r needs to be large to increase the potential D2D links. As K grows, a small 
r may result in a large number of the potential D2D links /cdb, which generates strong mutual 
interference to each other and limits the number of the scheduled D2D links. Thus, a smaller 
optimal r is needed for a larger K. Besides, we observe that the optimal r in our algorithm 
is smaller than the optimal r in [l26ll . This is reasonable since each user in our algorithm can 
establish more potential D2D links than the algorithm in ll26l for a given number of users. Then, 
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Figure 5. The Optimal file caching coefficient versus different file request coefficients, where dt = 0 dB, Cs = 0 dB, K = 100. 


a smaller optimal r is needed to maximize ks + ^D- 

Fig. Ogives the optimal 7 c with different file request eoeffieients 7 ^. For eomparison, we also 
provide the optimal 7 c in [|^ . It is observed that, the optimal r inereases as 7 c grows. When 
7 r is small, more different files are requested by the users in the eell. Then, users need to eaehe 
more different files to maximize the number of the seheduled D2D links. Thus, the optimal 7 c 
is small. Otherwise, a big optimal jc is required. Besides, we observe that the optimal 7 c in 
our algorithm is smaller than the optimal 7 c in [|2^ . This is beeause eaeh user in our algorithm 
ean potentially help more users than that in [[2^ . Then, users in our algorithm is required to 
eaehe more different files for D2D transmissions. Thus, a smaller optimal 7 c is needed in our 
algorithm. 

Fig. eompares the number of the seheduled D2D links of the proposed CPC-based D2D 
link seheduling algorithm with the DCPC-based D2D link scheduling algorithm in [[^ and the 
optimal schedule algorithm from the exhaustive search. Since the computational complexity of 
the exhaustive search with a large number of potential D2D links /cdb is too high, we limit the 
number of the potential D2D links by choosing a small help distance instead of the optimal 
one. It could be concluded from this figure that, the performance of the proposed algorithm 
is between the performance of the optimal schedule algorithm and the DCPC-based schedule 
algorithm in [[32l . This indicates that the CPC-based D2D link scheduling algorithm is more 
flexible than that the DCPC-based D2D link scheduling algorithm in [l32l . 

Fig. |7] shows the system throughput of the scheduled D2D links with the binary-search 
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Figure 6. Comparison of the numbers of the scheduled D2D links with different scheduling algorithms, where ut = 0 dB, 
Cs = 0 dB, 7 r- = 0.6, 7c = 1.5, r — j km. 



Figure 7. Comparison of the system throughput with different power allocation algorithms, where Cs = fx, K = 100, 7r = 0.6, 
7c = 1.5, and r = j km. 


based power alloeation algorithm in Algorithm |3l Here, the system throughput is ealeulated by 
log(r;m + !)• For comparison, we provide the performance of the optimal power allocation 

me<Sp 

algorithm, where a convex optimization toolbox is adopted, and the algorithm without optimiza¬ 
tion, where each D2D link works with the minimum SINK. From Fig. |71 the performance curve 
of the proposed binary-search based power allocation algorithm almost overlaps with the optimal 
one. This validates our analysis and indicates that our proposed algorithm may achieve similar 
rate performance with the optimal algorithm and outperforms the algorithm without optimized 
power allocation, i.e., more than 40% improvement. It should be noted that the proposed optimal 
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Figure 8. Performance comparison with scheduling algorithm in OH for different scheduling coefficient Cs, where ut = 0 dB, 
K — 100, 7c = 1.5 and r = 7 km. 


power allocation algorithm is based on max-min fairness, which achieves perfect fairness. In fact, 
the system throughput can be further improved if some unfairness can be tolerated If^ . 

Fig. [8] shows the number of the D2D links and the system throughput with different scheduling 
coefficient Cg. We observe that the number of the D2D links decreases as Cg grows. This coincides 
with the intuition that a larger Cg means a stricter condition on the scheduled D2D links. Then, 
the D2D links with stronger communication channels and weaker interference channels are 
scheduled. On the other hand, the system throughput increases as Cg grows. This is quite 
reasonable since fewer D2D links mean less mutual interference, which boosts the system 
throughput. On the other hand, the system throughput of the optimal scheduling in terms of the 
largest number of the scheduled D2D links is much smaller than the largest system throughput. 
Thus, the number of the scheduled D2D links and the system throughput can be balanced by 
choosing a proper Cg. 

Meanwhile, we provide the performance of the algorithm in lf34l . where only the information 
theoretic independent sets are scheduled. Note that there is no power allocation in f[34ll . we apply 
Algorithm |3] in [IM1 for fair comparison. From the figure, the number of the D2D links and 
the system throughput are constant since the scheduling algorithm in [f34l is not affected by Cg. 
We also observe that the number of the scheduled D2D links and the system throughput with 
the proposed algorithms can simultaneously outperform those with the algorithm in [l34l . For 
instance, Cg should be chosen around between four and six in this figure, e.g., Cg = 5 dB. In this 
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Figure 9. Performance comparison with the scheduling algorithm in 1341 in terms of fco and system throughput for different 
vt, where K = 100, jc = 1-5 and f = y km. 


way, we may choose proper values of Cg for different system parameters. 
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Table I 

The values of cs with different ut in Fig. [3 


Fig. |9] compares the proposed algorithm with the algorithms in ll34l in terms of the number 
of the scheduled D2D links and the system throughput. With the same selection method of Cg 
in Fig. [8l we choose the values of Cg with different vt in Tab. I. It is observed that the proposed 
algorithms outperform the algorithm in [f34ll in terms of both the number of the scheduled D2D 
links and the system throughput. 

Fig. [ini compares the proposed algorithm (Pro. in the figure) with algorithm A (Alg. A in the 
figure) and algorithm B (Alg. B in the figure) in terms of the number of the scheduled D2D links 
and system throughput. In this figure, (I) and (II) denote kjj and ks + kjj, respectively. In algorithm 
A, we adopt the Zipf-distribution caching and scheduling scheme in If26l . and the proposed 
optimal power allocation algorithm with perfect fairness meanwhile considering the minimum 
acceptable SINK constraint, i.e., ux = 0 dB. In algorithm B, we adopt the optimal cluster-based 
caching and scheduling scheme in [l25l . and the proposed power allocation algorithm with perfect 
fairness meanwhile considering the minimum acceptable SINR constraint, i.e., ux = 0 dB. The 
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K K 

Figure 10. Performance comparison with the scheduling algorithms in (251 and m in terms of ks -\- ko, kr>, and system 
throughput for different number of users in the cell, where it = 0 dB, ■jr = 0.6, 7c = 1.5, f = y km. 

scheduling coefficient Cs of the proposed algorithm is obtained with the same selection method 
in Fig. [8] for each K. That is, Cg is chosen to be 10 dB, 8 dB, 6 dB, 4 dB, and 2 dB when K is 
100, 150, 200, 250, and 300, respectively. From this figure, the proposed algorithms outperform 
the algorithms in Il25l and ll26ll in terms of the number of the scheduled D2D links, the number 
of users that can obtain the requested files either from their own memory or from helpers through 
D2D transmissions, and the system throughput. This is intuitive since the proposed algorithms 
create more D2D links with strong communication channels and weak interference channels. 
Then, the number of the scheduled D2D links and the system throughput can be simultaneously 
enhanced by efficiently scheduling and power allocation. 

Fig. [m compares the three algorithms (Alg. A, Alg. B, and Pro. in Fig. [TOl) in terms of the 
download time, which is calculated by EmeSBuSDB\S5 where is the 

download time of a video file from the BS to user m and is the download time of a video 
file from a helper to user m through an one-hop D2D transmission. Specifically, is calculated 
by where Tvideo is the length of a file, say Lvideo = 30 MB in the simulation, and 

is the transmission rate from the BS to user m and is assumed to be = 120 kbps 
[[261 , V m G S'b U S'db\ 5'6. Similarly, is calculated by = 4^, where R^^ is the 
transmission rate of D2D link and is assumed to be R^^ = log(um + 1) Mbps. From this 
figure, the proposed algorithms outperforms both Alg. A and Alg. B. This is reasonable since 
more users in the proposed algorithms can obtain the requested files from their own memory or 
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K 

Figure 11. Performance comparison with the scheduling algorithm in l26l in terms of the download time for different number 
of users in the cell, where ut = 0 dB, 7 ^ = 0.6, 7 c = 1.5, f = 7 km. 


from helpers through D2D transmissions eompared with those in either Alg. A or Alg. B (as 
shown in Fig. [TOl). Thus, the average download time is redueed with the proposed algorithms. 

VII. Conclusions and Future Discussions 

We have studied the effieient seheduling and power alloeation of D2D-assisted wireless eaehing 
networks. We formulate a joint D2D link seheduling and power alloeation problem to maximize 
the system throughput. However, the problem is non-eonvex and obtaining the optimal solution 
is eomputationally hard. Alternatively, we seek to obtain a suboptimal solution with reasonable 
eomplexity. Briefly, we intend to sehedule the D2D links with strong eommunieation ehannels 
and weak interferenee ehannels and alloeate the power to the seheduled D2D links fairly. Thus, 
we deeompose the system throughput maximization problem into a D2D link seheduling problem 
and an optimal power alloeation problem. To solve the two subproblems, we first develop a D2D 
link seheduling algorithm to inerease the number of the seheduled D2D links satisfying both 
the SINK and the transmit power eonstraints. Then, we develop an optimal power alloeation 
algorithm to maximize the minimum transmission rate of the seheduled D2D links. Numerieal 
results indieate that both the number of the seheduled D2D links and the system throughput ean 
be improved simultaneously with the Zipf-distribution eaehing seheme, the proposed D2D link 
seheduling algorithm, and the proposed optimal power alloeation algorithm eompared with the 
state of arts. 
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From the results in this paper, we also conclude that the D2D communication in a wireless 
caching network can be enhanced by creating more D2D links with strong communication 
channels and weak interference channels. This can be achieved by optimized caching, the 
proposed D2D links scheduling algorithm, and the proposed optimal power allocation algorithm. 
Although the optimized Zipf-distribution caching scheme in this paper is not generally optimal, 
the proposed D2D link scheduling algorithm and the proposed optimal power allocation algorithm 
can be used for any caching scheme. Ideally, the joint caching, scheduling, and power allocation 
algorithm should be able to further improve D2D communication. However, even the optimal 
caching distribution of the general model in our paper is quite hard. A further study on the 
optimal caching scheme is beyond the scope of this paper. 

Throughout this paper, we assume one channel to simplify the illustration of our algorithms. 
In fact, the proposed D2D link scheduling algorithm and the proposed optimal power allocation 
algorithm can be used in multi-channel systems, e.g., FDMA and TDMA. More specifically, 
D2D links are iteratively scheduled in a unique channel, i.e., a frequency band (corresponding 
to FDMA) or a time slot (corresponding to TDMA). This process terminates in two cases. One 
is that there is no more channels to accommodate the scheduled D2D links. The other one is that 
all the D2D links have been scheduled. Since more wireless resource is used in multi-channel 
systems, a better system performance is expected. 

Besides, both the proposed D2D link scheduling algorithm and the proposed optimal power 
allocation algorithm are achieved in a centralized manner. That is, all the calculations are 
conducted at the BS. This requires the BS to keep track of the cached content in each memory 
and collect the channel state information among different potential D2D links. One direction of 
the future work is to develop decentralized algorithms of scheduling and power allocation in the 
D2D assisted wireless caching networks. 

VIII. Appendix 

A. Proof of Theorem 1 

To prove Theorem I, we will first prove part I: each transmit power pn,^ n e S^, increases if 
any SINK v^, rn G iSd, in (Pi) increases. Then, we will prove part II: the value of |iSd| remains 
constant or decreases if any SINK Vm, in G iSd, in (Pi) increases. 



31 


1) Proof of part I in Theorem 1: Suppose that there is a power inerement and a SINK 
increment for D2D link in iSd such that 


Vm + Ap 


{Pm + Apm)g{m,m) 

EneSo,n^mPn9{n,m) + N^' 


( 22 ) 


Then, there must be a SINK decrement Avn for any other D2D link In (Vn E SY),n m) such 


that 


Vn - An 


n 


_ Png(n,n) _ 

J2k&Su,k^n,k^mPk9(k,n) + (Pm + Ap^)c/(m, n) + 


(23) 


To remain the SINK n„ at D2D link In, Pn ('^ n E iSd, n m) will be increased by Ap„, such 


that 


_ (Vn + Apn)g(n,n) _ 

J2k(,Sn,kMPk + ^Pk)9{k,n) + Nn 


(24) 


Thus, each p„, n G iSd increases if there is an increment of any Vm, m E Sj). 

2) Proof of part 11 in Theorem 1: Suppose that the D2D links in £d satisfy all the constraints 
in problem (Pi) with transmit power P = {pi,p 2 , • • • ,Pko}- If there is a SINK increment Avm at 
D2D link 1^, the transmit power should be increased to P+AP = {pi+Api, p 2 +Ap 2 , • • • , Pko + 
Apkjj} to satisfy all the minimum acceptable SINRs. If pn + Apn < 'i n E iSd, all the 
D2D links in £d can still be satisfied with the minimum SINK constraints and number of the 
scheduled D2D links remains. If there exists a Apn (n E S^) satisfying pn + Apn > pf‘^^ and 


pf^^ 9 (n, n) 

Efce5D,fc^n(Pfc + ^Pk)9{k, n) + Nn 


(25) 


Then, the D2D link In cannot be satisfied with the minimum SINK constraint. This reduces the 
number of the scheduled D2D transmissions. This completes the proof of Theorem 1. 


B. Proof of Theorem 2 

The regular admission control problem with QoS constraint has been proved to be NP-hard 
in [|3^ and is a special case (Case II) of problem (P 2 ). Specifically, the problem (P 2 ) reduces 
to a regular admission control problem when any two D2D links in do not share the same 
users. Thus, problem (P 2 ) is also NP-hard. 
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C. Proof of Theorem 3 


In this part, we will prove Theorem 3. Firstly, we will give the proof of the first ease that 
there is a user m (0 < m < /cd) satisfying P ^ ^ but not satisfying 

P ^ P(V("*)(r;m+i)) ^ pmax^ prove that the optimal SINRs of the seheduled D2D links 
is V* = [ u*, • — , u* , Vm+i, ■ ■ ■ where Vm < v* < Vm+i- Then, we will provide the proof 

m _ 

of the seeond seeond ease that if P ^ ^ pmax bolds, the optimal SINRs of the 

seheduled D2D links is V* = \ v*, ■ — , v*j , where V* > Vk^. 

kj) 

1) Proof of the First Case: We will prove this by two steps by eontradietion, the first step is 
to prove the first m optimal SINRs are identieal, i.e., u* = u* for 1 < n < m. The other step 
is to prove the last k^ — m optimal SINRs are the minimum SINR eonstraints, i.e., u* = Vn for 

m + 1 < n < /cd- 

Denote the optimal SINRs of the seheduled D2D links as V* = [ u*, • ' ‘ ‘ 

Sinee P ^ ^ p™^ holds and P ^ P(V*^™')(um+i)) ^ p™^ does not hold, we 

have Vm < < Vm+i for 1 < n < m and u* > for m + 1 < n < fco- If we denote the 

optimal power alloeation as P* = [pfp 2 , ■ ■ ■ ,pIJ, we have 


Pn9{n,n) 


Ellrk^nPl9{k,n) + N„ 


= > Vm > Vn,y I < n < m. 




(26) 


We observe that the objeetive funetion in (P 3 ) is equivalent to min{ max — | and suppose that 
the optimal power alloeation results in different SINRs at the first m D2D links, i.e., max 1/u* > 

l<n<m 


min 1/v* and that the D2D link In* has the largest SINR at first m D2D links, i.e., n* = 

l<n<m 

arg min 1 /u*, wehave vl* > Vn,. Besides, from (l26l) . we observe that u* is a strietly inereasing 

l<n<m 

funetion of Pn and is a strietly deereasing funetion of pi for k n. Therefore, there must be 
a small power deerement and SINR deerement for D2D link In* and a small SINR 
inerement Au„ for other D2D links In { ^ ^ n < ku, n n*) sueh that the eonstraints in (P 2 ) 
still hold, i.e.. 


{p*n, -/\p)g{nfn*) 


E n=kD 
n=l,n^n 


.p*ng{n,n*) + Nn* 


= 


- ^Vn* 


A Vni ^ Vn* 


(27) 


and 


p*ng{n,n) 


Y.sli,s^n*,s^sP*s9{s,n) + {pl 


/\pn*)g{n*,n) + Nk 


= V*n + AVr, 


> Vr, 


(28) 





33 


Suppose that we choose a small to satisfy 
have 


1 

(<* ) 


< max 

l<n<m, n^n* 


1 

{v*+Av„) 


. Then, we 


max 

l<n<m 



max < 

rii 

> > max < 

l<n<m, n^n* 


l<n<m, n^n* \ 




(29) 


= max < max <-r— > ,-r-> . 

(l<n<m, n^n* ( W* + J V** — Av^* J 

This means that there exists another power allocation P' ^ P* enabling the SINRs at the 
scheduled D2D links to be V' ^ V* and max I/7* > max 1/v'^, where 1/v'^ is equal to 

l<n<m l<n<m 

max + Avn) for n ^ n* and is equal to — An„*) for n = rf, which causes 

l<n<m, n^n* 

contradiction. Thus, we have max l/n* = 1/v*, where Vm A v* < Vm+i- 

l<n<m 

Then, the optimal SINRs vector can be denoted as V* = [n*, • • • and the 


corresponding minimum SINR is v*. Next, we will prove n* = for m + 1 < n < ku- 

Suppose that there exists a user k* satisfying for m + 1 < /c < ku- There must be 

a small power decrement Apk* and SINR decrement Avk* for D2D link Ik* and a small SINR 
increment Avk for D2D links h ( ^ A k < m) such that the constraints in (P 2 ) still hold, i.e.. 


{P*k* - Apk*)g{k*,k*) 


E t=ki) 

k=l,k^k 


*P*k9{k,k*) + Nk* 


= vl. - Avk* > Vk*, 


(30) 


and 


_ P*k9{k,k) _ 

T.t=i,t^k,t^k* P*t9{t, k) + {pl, - Apk*)9{t, k) + Nk 


V* + Avk- 


(31) 


Then, if we choose a small Apk* and Avk* to satisfy min v* + Avk < Vm+i, the minimum 

l<k<m 

SINR of the scheduled D2D links is min v* + Avk, which causes contradiction. Thus, we have 

l<k<m 

n* = Vn for m + 1 < n < /cd- 

This complete the proof of the first case. 

2 ) Proof of the Second Case: The proof of the second case is similar to that of the first case 
and will be omitted for page limit. 
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